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Abstract. We prove a maximum principle for the problem of optimal control 
for a fractional diffusion with infinite horizon. Further, we show existence of 
fractional backward stochastic differential equations on infinite horizon. Fi- 
nally, we illustrate our findings with an example. 



1. Introduction 

In this paper we consider a control problem with respect to a given performance 
functional: 



J(u) = E 



f(t,X(t),u(t),w)dt 



o 



where X(t) is a controlled fractional diffusion and u{t) is the control process. We 
allow for the case when the controller only has access to the partial information 
£ t at time t. Thus, we have a infinite horizon problem with partial information. 
Infinite-horizon optimal control problems arise in many fields of economics, in par- 
ticular in models of economic growth. Note that because of the general nature of the 
partial information filtration £ t , we cannot use stochastic dynamic programming 
and the Hamilton- Jacobi-Bellman (HJB) equation to solve the above optimization 
problem. Thus our problem is different from partial observation control problems. 

In the deterministic case the maximum principle by Pontryagin (1962) has been 
extended to infinite-horizon problems, but transversality conditions have not been 
taken into account in the litrature in general. The 'natural' transversality condi- 
tion in the infinite case would be a zero limit condition, meaning in the economic 
sense that one more unit of good at the limit gives no additional value. But this 
property is not necessarily verified. In fact |Hal74| provides a counterexample for a 
'natural' extension of the finite-horizon transversality conditions. Thus some care is 
needed in the infinite horizon case. This paper contains a extension of the maximum 
principle in [BH0SO2] and the existence results of fractional backward stochastic 
differential equations found in |HP09j to infinite horizon. 



In the case of Brownian motion there have been several paper on infinite-horizon, 
see e.g. |H0P12| , and |HM0P12| and (AH0P12| for the delay case. We are going 
to extend the results in [H0P12J to the case of a controlled diffusion process driven 
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by a noise with memory. Here the driving noise in our controlled system with 
memory is modelled by a fractional Brownian motion, B^ 1 , t > 0, where the Hurst 
parameter h € (i, 1). By solving this problem we establish a sufficient stochastic 
maximum principle. The latter result requires the construction of unique solutions 
of fractional noide driven backward stochastic differential equations (fBSDE's) of 
the form 

jdp(t) = -f(t,r) t ,y t ,z t )dt- z t dB^ , 
] limp(i) = 0. 

V t— ► oo 

As for the theory of BSDE's, which has gained more and more importance in 
stochastic control and mathematical finance, we refer the reader to the seminal 
paper by Pardoux and Peng |PP90| . For some technical background to fractional 
Brownian motion, see |Ben03| and [BH0ZO8|. The latter one will be a basis for 
much of our framework. For a comprehensive background to the mathematical 
structures used, see [BK95J. 

2. Preliminaries 

Let B" , t > 0, be a fractional Brownian motion with Hurst index H £ (i, 1) on 
the probability space (fi, F H , P H ). We assume (fl, F H , P H ) is endowed with the 
natural filtration of B H , where F H = VtyoJ-^ 1 ■ 

In the following we aim at introducing stochastic integrals with repsect to B^ 1 
by using techniques from Gaussian white noise analysis. For a general introduction 
to White noise and Hida distributions in the Brownian motion case, the reader 
may consult the excellent books [H0U Z1O] , |Kuo96| . |Oba94| and of course the 
classical book [HHPS93J. For a comprehensive explanation of the contruction of 
the fractional Brownian motion, where many of the following theorems and lemmas 
are found, we refer to )BH0ZO8| . 

Let 1 > H > i, and 

t) = 4> H {s, t) := H(2H - t\ 2H - 2 . 

Denote by <S(R) the Schwartz space of rapidly decreasing smooth functions on R, 
and for f,g £ S(M), define 

(f,9)H,t-= / f{s)g{r)(t){s,r)dsdr, 
Jo Jo 

and 

(/,<?) H := / f f(s)g(t)<f>(s,t)dsdt. 



Now consider the completion of S(R) under this inner product and denote the 
resulting separable Hilbert space by 

(1) L% 

Let n — S'(M) be the topological dual of <S(R) (the space of tempered distributions). 
By Bochner-Minlos theorem or a abstract Hilbert space argument we have that 
there exists a probability measure P H on the Borel cr-algebra such that 



e 



(»>f)dF H (cj) = e-iM», 
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for all / € <S(R), where (•, •) is the dual pairing between <S(R) and <S'(R). It is now 
easy to see from a Fourier argument that 

Epb [<•,/)] :=£[<•, /)]=0, 

and 

E[{;f) 2 ] = \\f\\ 2 H . 

Now we define the fractional Brownian motion as 

Bf = B H (t,u) = {u,l [0>t] (-)}, 

by extension, which is an element of L 2 (V H ) for each t. Let J = (N{J) C be the set of 
all finite sequences a — (ai, . . . , a m ) of nonnegative integers. For a = (ai, . . . , a m ), 
let 

H a (uj) := h ai ((u},ei))---h am ((u,e m )), 

where h n is the n-th Hermite polynomial (see [H0UZ1O] ) and {e n } n >o a orthonor- 
mal basis of L^(W). As in the case of Brownian motion we obtain a Wiener- ltd 
chaos expansion. 

Theorem 2.1. let F e L 2 (P H ). Then there exist unique C a , ot € J such that 

F(u) = c a H a (uj), 

where the convergence is in F e L 2 (P H ). We also have that 

\\F\\ 2 h= E a!c - 

where a! := ax ■ ■ ■ a m . 

See IH0 UZ1Q] , Theorem 3.1.8 for a proof. We are now ready to define the 
fractional Hida test function and distribution spaces. 

Definition 1 (The fractional Hida test function space). Let (S)h be the set of all 

= E aeJ «^«H e L 2 (F H ) such that 

where 

(2N) Q =n(2j) 73 - 

3 

Definition 2 (The fractional Hida distribution space). Let (S)* H be the set of all 
formal expansions 

G(u) = E b p Hp(u) e L 2 (P H ) 

such that 

\G\\k- q = E wimr q p, 

for some q E N. 
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Let (S)h be equiped with the projective topology and (S)* H with the inductive 
topology. Then (S)* H is isomorph (in the category of topological vector spaces) to 
the topological dual of (S)h with action given by 

<€.G,tp >:= (G, ip)( S ) H ■= a - a ^- 

Definition 3. Let 

F(oj) = ^ a a l-L a {uj) 

and 

G(w) = hKp{Lo) 
belong to (S) * H . Then we define the Wick product F o G by 

FoG(w) = dabpUa+p^) 

a,/3ej 

Further, let 

DfF= f D?F<t>(t,s)ds 
denote the Malliavin (^-derivative of F (see [BH0ZO8J Definition 3.5.1). 

Definition 4 (Fractional Wick-lto-Skorohod integral (fWlS)). Suppose Y : R — >• 
(S)g is a function such that Y(t) o W H it) is Bochner (or weaker Pettis) integrable 
in (S)* H . Then we define the fWIS-integral, J R Y(t)dB^, as 

/ Y(t)dB? := / Y(t) o W t H dt 
Jr Jr 

3. Fractional Wick-Ito-Skorohod integral in L 2 

The integral defined above is an stochastic distribution, but we would like to 
work in L 2 so we need to do some work to ensure that. To achieve this, we extend 
the construction of the fractional Brownian integral in [BH0ZO8J to infinite horizon. 

3.1. Integral of simple functions. Consider a partition 7r„ of [o,n], it : = to < 
t\ < . . . < t n = n. We now define the integral 



i 

for a a simple function F n £ L 2 (F H ) of the form 

n 

r(t,,) = ^^Hi [Wi+l) (t), 

as 

/ F^dBf := S(F, tt„) := £ F u o - B*), 

if |j/|| £ i,2 < oo where 



|/U, := -B 



f(s)f(t)0(s,t)dsdt + ( Dff(t)dt) 
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Using that for F, G e L 2 (P H ) we have that E[F o G] = E[F]E[G] and 

n r 

E[S(F,* n )\ = Fr(u)»(B* +i -B*) 



i=0 



= J2e ftwe{b»^-b%) 



»=o 



= 



Further we have that 



E 



(FP{u) o {B« +i - B«))(F?{u>) o (B« +1 - B%)) 



D*F£BfF£dtds + F^F i 



4>{s, t)dsdt 



so that 



E[S(F n 1 7r n ) 2 ] = E 

i,j=0 

+ *K 



rti+i rtj+i 
Jti Jti 



tj+i rtj+i 



D f F n D f F n dtds 



U 



<j)(s, t)dsdt 



3.2. FWIS integral for general stochastic functions. Define 

\ir n \ := max(t i+ i - U) 

i 

and := / tj if U < t < t i+1 . Assume that E[\\f^ - f\\ 2 H ] -t and 
E[S(F, 7r n ) 2 ] =J2 E [ + [ DfF u DfF u dtds 

i,j=0 *-Jti Jtj 



+ F U F U 



ti+l ftj + 1 



<fi(s, t)dsdt 



converges to as |7r n | — > 0. For a sequence of partitions {7r„} such that |7r n | — > 
then S(F,n n ) is a Cauchy sequence in L 2 (¥ H ). The limit in L 2 (F H ) is 



(2) 

so that 



/ /(*)<iB* (5) := lim ^ / ti o (B* +1 - B*), 



/ /(*)dB"( fi 

Jo 



2 








= £ 








. Jo Jo 



DtF u DtF u dtd S + \\f\\ 2 H 



Definition 5 (£^,(0, oo) ). Let £^(0,00) be the family of stochastic prosesses, /, 
on [0, 00) with the following properties: 

(1) E\\f\\% < 00, 

(2) / is 0-differentiable, 

(3) the trace of Dff t , < s < t < 00 exists, 

(4) E[f °°f °°\Dtf t \ 2 dsdt}<™ 
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(5) and for each sequence of partitions (ir n ,n € N) such that |7r„| — > as 
7i — y oo we have that 

*U+i rtj+i 



i,j=0 

and 



£>[ / / 3+1 \D*F§D+F& - D*FZD+FZ\dtd8 

„_n Jti Jti 



E[\\F* - F\\ 2 H ] 



tend to as n — > oo, where n n := = < t" < . . . < t'' l \_ 1 <t™ = n. 

We summarize the above in a theorem. 

Theorem 3.1 (Properties of the integral). Let {ft}t>o be a stochastic process such 
that f 6 £,7,(0,oo). Tfte Hmii (H]) exists and satisfies; 

/>oo poo />oo 

E[( f{s)dB H {s)) 2 ]=E[ / D*( a )/(t)D*(t)/( a )d*ft+||/||ff] 
Jo Jo Jo 

and 

/•OO 

E[ f(s)dB H (s)}=0. 
Jo 

3.3. Itos formula. In this section we present Itos formula for fractional Brownian 
motion. Fist let us look at the </S — derivative: 

Theorem 3.2. Let (Ft,t € [0, 00)) be a stochastic process in £^([0,00)) and 
sup E[\DfF s \] < 00 and let r} t = t G u du + f* F u dB" . Then for s,t>0 

0<s<cx> 

ft 

F u 4>(s,u)du, 



D%= [ DtG u du+ f DfF u dB%+ [ 
Jo Jo Jo 



For a proof see [BH0ZO8J. Now let we arrive at Itos formula: 

Theorem 3.3 (ltd formula). Let rj t = 770 + Jq G u du + f F u dB^ , where (F t ,t € 
[0,oo)) is a stochastic process in £^([0, 00)). Assume there is an a > 1 — H such 
that 

E[\F u -F v \ 2 ]<C\u-v\ 2a , 
where \u — v\ < S for some S > and 

Jim E[\D£(F u -F v \ 2 ] = 0. 

0<«,u<t,|«— v\— >0 

Let f : R + xR — > R be a function with continuous derivative in the first variable and 
twice differ entiable with continuous first and second derivatives in the second vari- 
able. Assume all derivatives are bounded. Moreover, assume E[J^° \F s Dfr) s \ds] < 
00 and (/ (s,ri s )F s ) is in £^([0, 00). Then 

f(t, Vt ) =/(0,0)+ / ^f(s, Vs )ds+ I ?f(s,r) s )G s ds+ [ ?f(s,r) s )F s dB? 



OS Jq OX Jo ox 

1 8 2 f 

J (s,r] s )F s DU s ds 



lo dx 2 
s. 

For a proof see [BH0ZO8] . 
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3.4. Fractional Clark-Hausmann-Ocone theorem. In this section we briefly 
recall a Clark-Hausman-Cone theorem for fractional Brownian motion on a distri- 
bution space L 2 (P H ) C Q* C (S)* , which we want to employ in Section 4. Let us 
first give a construction of Q* to which the operator D t will be extended. 



Definition 6 f |PT95| . [A0PUOO|4BH0ZO8| ). Let k e N . We say that a 
random variable with chaos expansion 



^ = £ / fnd(B H r n (t), 

.,_n JR*! 



,i=0 " M + 

where /„ e i^(R™ ) belongs to the space G k = G k {P H ) if 



k ■.= J2 nl \\f"\\Ll(M V e 2kn <oc, 
where ||-||l2( R ,ii is the completion of 5((R") as in (JXJ) . Now. letting 



g = G(r H )= n g k (F H ) 

k>0 

and equip Q with the projective topology. 
Let q 6 No- We say that a function 



oo „ 



n=0 

where # n € L|j(R!j:) belongs to the space Q- q = Q_ q {F H ) if 

where ||/n||z, 2 (R«) is the completion of <S((R™ ) as in ([T]). Now, letting 

g* =g*(p H ) = u g- q {P H ), 

q>0 

and equip g with the inductive topology. Then g* is the dual of g, and the 
action oi G £ Q* on ip £ Q is given by 



rWL 2 ,(R")- 



n=0 



We will also need a variant of the conditional expectation on g* , that is easier 
to work with. 

Definition 7. (1) Let 

oo „ 

G = J2 9n(s)d(B H r n ( s )eg*. 

Then we define the fractional conditional expectation of G with respect to 
Tf by 



£[F|jf ] = / . 9 „(,s)l {0 < s < f} d(B ff )^( s ) 



for t > 0. 
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(2) We say that G € Q* is J" t ff -measurable if 

E[G\T?] = G, 

for t > 0. 

Remark 3.1. The quasi-conditional expectation E^J 7 ^] was introduced in [H0O3 
to construct hedging strategies of financial claims. 

An immediate consequence of the definition is the following, which is important 
in solving backward stochastic differential equations. 

Lemma 3.4. If (f u ,0 <u<T) is a real valued stochastic process such that 

rp rp pi p 

E[ [ <t>(s-t)\f t \\f s \dtds+ [ Dfftdsdt], 
Jo Jo Jo Jo 

and jf fudB? 6 L 2 (P H ) then 

E[( f u dB*\F t ] = 0. 
Jo 

Remark 3.2. Note that this also holds for T = oo. 

Proof. Let f u = X^^Lo In{gn(u)) In be the chaos expansion of / , where g is a 
function of n- variables. Then it is well known that 

(3) / fudB% = (sn(u)), 

where g n +i is the symmetrization. Now the result follows from the definition of the 
quasi-conditional expectation. □ 

Lemma 3.5 (Clark-Hausmann-Ocone representation). Let F € L 2 (P H ) then 

/•OO 

F = E[F] + / E[D t F\T t H ]dB? 
Jo 



Proof. 



OO r. 

= E[F] +J2 ■ • • » »»)dB£ • • • dB£ 

n= l *^[0<si,--- ,s T1 <oo] 

°° /"OO n 

= E[F] + Y / (/ f(si,...,s n - l ,t)dB? 1 ---dB? n _ i )dB* 

n=1 Jo J[0<s u -- ,s n -i<t] 

°° /"OO 

= ^ + E / ^* 7 " (/«) \?t]dB? 

= E[F] + / E{D t F\T t H ]dB?. 
Jo 



□ 
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4. Fractional backward stochastic differential equations (fBSDE) 

Let -Bj , t > 0, be a fractional Brownian motion with Hurst index H > h on 
the probability space (f2, F, P ff ) endowed with the natural filtration of and 
T = U t > J"f. Let 6 : [0, oo) xlxl — > R. Consider the problem of finding a 
J^-adapted processes (p(t),q(t)) such that 

Jdp(t) = b{t,p(t),q(t))dt + q(t)dB H {t), 

| limp(t) = 0. 

V i— Voo 

This is a infinite horizon fractional backward stochastic differential equation (ihfB- 
SDE). 

4.1. Existence of general FBSDE. In this section we prove a result about exis- 
tence and uniqueness of the solution (Y(t), Z(t),K(t, £)) of infinite horizon BSDEs 
of the form; 

(4) dY(t) = -g(t,Y(t),Z(t),K(t,-))dt + Z{t)dB H (t);0 < t < r, 

(5) limy(t)=C(r)l [ o, 00 )(r), 

where t < oo is a given J-^-stopping time, possibly infinite. Our result is the 
fractional vesion of [H0P12J, and the infinite horizon version of |HP09j . See also 
|Par99| . |Pen90| . |Yin08| . [LT09] . [LT94] . |Roy06| , |BBP09j and [Slt02j . for the 
classical Brownian motion case. 

Let rjt = r/o + b t + / a s dB^ , where 770 is a given constant, 6 t is a deterministic 
difierentiable function of t, and cr s is a deterministic continuous function such that 
<j t exists for all t and ^||c||t exists and it is strictly positive. Let £ = h(ri T ) and a 
continuous h . Let 

1 _*£ 
V Z7TI 

Let the semigroup {Pt}t>o be given by 

f-*/(aO = I Pt(x-y)f(y)dy. 

JR 



If f is continous we have that 



f( Vt ) = P M 2f(m) + I £iVnM k || ? /fa>(a)dBf • 



Using Ito, we get 
Therefore 
Define 

Voo = {Y(-) = 4>(; 77(0; G C^ 2 for all t G [0, 00)}, 

where C { ' is the set of functions that are continuously difFerentiable with respect to 
t and twice continuously difFerentiable with respect to x. Let Voo be the completion 
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of Voo under the norm 



/•OO />oo 

\Y\\*:= e xt E[\Y t \ 2 }dt = e xt EMt, Vt )\ 2 ]dt. 
Jo Jo 



We require the folowing: 
(1) For each t > T > 



inf — > Co, 

0<s<Ta* 



for some posetive constant Co, where 

6(s, r)a r dr. 



(2) The function g : ft x R + x 1 x K -> R is such that there exist real numbers 
fi, A and K, such that K > and 



2 



(6) A > 2^ + 2^ 

We also assume that the function g satisfies the following requirements: 
(a) g(-, y, z) is progessively measurable for all y £ R, z £ K, and 

\g(t, y, z) - g(t, y, z')\ < K\z - z'\. 

(7) 

(b) 

(8) (y - y', g(t, y, z, k) - g(t, y', z, k)) < n\y - y'\ 2 

for all y, y' , z a.s. 



(9) 



(c) 

E [ e xt \g(t,Q,0)\ 2 dt < oo. 
Jo 



(d) Finaly we require that 

(10) y^g(t,y,z), 

is continuous for all t, z a.s. 
(3) Let r)t = t)q + b t + J Q a s dB^ as above. Assume the final condition £ is given 
by £ = r] T , and f t = B[£|Ji] such that E(e Xr |f| 2 ) < oo and 

(11) £ / e xt \g(t,£t,(Tt)\ 2 dt <oc. 

Jo 

A solution of the BSDE ((4])-([5|), is a dual (Y t , Z t ) of progressively measurable pro- 
cesses with values in 1 x i s.t. Z t = when t > r, 

(1) £[/ r e As £>f Z s ds] < oo, 

(2) Y t = F TAT + J tA ^ T g s ds - J tA ^ T Z s dB^ for all deterministic T < oo and 

(3) Y t = £ on the set {t > r}. 

Remark 4.1 (Infinite Horizon). This incorperates the case where t(lu) — oo on some 
set A with P(A) > 0, possibly P(A) = 1. 
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Theorem 4.1 (Existence and uniqueness). Under the above conditions there exists 
a unique solution (Yt,Zt) of the BSDE ((4|)-([5|), which satisfies the condition; 



(12) 



E[[ T e Xs (\Y s \ 2 + \Z s \ 2 ds] 
Jo 

<cE[e XT \e+ [ T e Xs \g(s,0,0)\ 2 ds], 
Jo 



for some positive number c. 

Proof. First, let us show uniqueness: 

Let (Y,Z) and (Y',Z') be two solutions satisfying $TZ§ and let (Y,Z) = (Y - 
Y',Z — Z'). From Ito's Lemma we have that 



e AMT |^Ar| 2 
TAt 



< - 



e Xa \\Y a \ z ds + e AS \Y T \ z + 2 



As 



tAT 

TAt 



TAt 



Y s f(t,Y s ,Z s )ds 



e Xs DfY s Z s ds - 2 



tAT 



TAt 



e Xs Y s Z s dB s 



H 



<e As |F T | 

cTAt 

+ 



2\\Y S \ 2 + ^\Y S \ 2 + K\Y S \\Z S 



ds + 2 



TAt 



e Xs DfY s Z s ds 



TAt 



e xs Y s Z s dB s 



H 



tAT 



< e Xs \Y T \ 2 



TAt 



.As 



2\\Y S \ 2 + n\Y s \ 2 + K\Y S \\Z S \ 



ds — 2cq 



TAt 



e Xs \Z s \ 2 ds 



tAT 



fTAT 



e As Y s Z s dB s 



H 



The last inequality follows for the fact that for Y t = (j>(t,rit) we have that Z t 
—&t<f>x(t,i]t) (see |HP09| Proposition 4.3) so that 



DfY s = [ 0(s-r)D r Y s dr 
Jo 

= 4> x (s,ri s ) 4>(s-r)a r dr 



So 



= o^ s (t) x (s, r/ s ) 

= — -z s . 

OS 



DfY s = -—Z t . 
o~t 
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Then it follows that 



e XtAT \Y tAT \ 2 



TAt 



e Xs 2 



r TAT 

<e x < TAT >|Y T | a + / 2 

J t At 



A|r s | 2 + Co |z s | 2 js 

e Xs v\Y s \ 2 ds + K\Y s \\Z s \) 



ds 



TAt 



e Xs %ZAB". 



Combining the above with the fact that 2ab < jt- + 2c b 2 and A > 2/i + 2K 2 , 
we deduce, that for t < T 



e Xt ^\Y tAT \ 2 

/•TAt 

< e ( TAr )|r T | 2 -2 / 

JtAT 



e""YZAB". 



Letting T — > oo, so that we have, since E[J e At |F(i)| 2 ] < oo 



e XtAT \Y t A T \ 2 



< -2 I e As Y s Z s dB s 



H 



Taking expectation the uniqueness follows. 
Proof of existence: 

For each n G N we construct a solution (F t n , Z™) of the BSDE 

/>nAr t>T 

Y t n = £ + / g(s, y s \ z:)d s - / zyB B 

JtAT JtAT 

by letting {(Y t n , Z r t l );0<t< n} be defined as a solution of the following BSDE: 

l l0 , T] (s)g(s,Y s n ,Z2)ds- J Z?dB s 
for < t < n and {(Y t n , Z?);t>n} defined by 



and 



for t > n. From Theorem 4.6 in [HP09] we we now that this finite horizon equation 
has a unique solution thanks to our requirements. Next, we find some a priori 
estimates for the sequence (Y n ,Z n ,K n ). For any e > and p < 1 we have for all 
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t > 0,y e R, z e K with c=\, 

2(y, g{t, y, z)) = 2{y, g(t, y, z) - g(t, 0, z)) 
+ 2(y,g(t, 0, z) - g(t, 0, 0)) + 2(y, g(t, 0, 0)) 

< (2p+- p K 2 + e)\y\ 2 + p\z s \ 2 

+ c\g(t,0,0)\ 2 . 
From Ito's Lemma and arguments given above, we have 



e AtAT| y n r |2 + 



|2 + p|2 s | 2 l ds 



f T e As [\\Y?\ 

JtAT 

< e Xs \ri\ 2 + c [ e Xs \g(s,0,0,0)\ 2 ds 

JtAT 

-2 [ T e Xs <Y s n ,Z^dB s >, 

JtAT 

with A = A - 2(i — - p K 2 — e > and p = 1 — p > . From this it follows that 



E t 

I sAt 

< CE 



e Ar (|y r T+ || z? || 2 ) 



dr 



M\ 2 + I' e Xr \g(r,0,0)\ 2 dr . 

J sAt 

Let m > n and define AY t := Y t m -Y t n and AZ t := Z™-Z?, so that for n < t < 

pmAr ptuAt 

AY t = g( S ,Y s m ,ZT)ds- AZ s dB s . 

JtAT JtAr 



It then follows that 

cmAT 



/•THAT 

/ {e As (A|AF s | 2 + |AZ 5 | 2 ) 

JtAT L 

pmAT 

/ {e As (A|Ar s | 2 + |AZ s | 2 ) 



< e AtAr| A y tAr |2 + 



-X 



/■mAr 
JtA 



' t At 



/•mAr pmAT 

2<e Xs \j 1 \ 2 c e A %(s,0,0)| 2 ds-2 / e As Ay s AZ,dB s . 

JtAT JtAT 



From the same arguments as above 



E 



pmAT 

/ r 

JnAT *- 



I AY. 



< 4£ 



/ e As | 5 ( S ,^a)| 2 rf S . 

.JnAT 

The last term in the above equation goes to zero as n — > oo. Now, for t < n 

pnAT pnAT 

AY t = AY„ + / {«,(«, Y™, ZD - IT, Z?) U - / AZ s d£ s . 

JtAT L J JtAT 
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Using the same argument as in the case of uniqueness, we have that 



£[e AtAT |Ay tAT | 2 ] < E[e MAr \AY n \ 2 } < cE 



e As | ff ( S ,£ s ,77 s )| 2 ds 



It now follows that the sequence (Y n , Z n ) is Cauchy in the norm 

\\{Y,Z)\\ e Xs {\Y s \ 2 + \Z s \ 2 )ds]. 

Jo 

So, we have that there is an unique solution to the BSDE ((4j)- ((5j) , which satisfies 
for all A > 2/i + 2K 2 , the condition 



e\ T e Xs (\Y s \ 2 + \Z s \ 2 )ds 
Jo 



< cE 



Jo 



e Xs \g(s,0,0)\ 2 ds 



□ 



4.2. linear infinite horizon backward stochastic differential equations (LI- 

HBSDE). For linear infinite horizon backward stochastic differential equations 
(LIHBSDE), we can give a constuctive proof of existence. 

dp(t) = [a(t) + b(t)p(t) + c{t)q{t)]dt + q(t)dB H (t), 
lim p(t) = 0. 



(13) 



where b(t) and c(t) are given continuous deterministic functions and a(t) = a(t,u>) 
is a given F^-adapted process such that J \a(t,u)\dt < oo almost surely. By the 
fractional Girsanov theorem we can rewrite (|13D as 

(14) dp(t) = [a(t) + b(t)p(t)]dt + q{t)dB H (t), 

where 

B H {t) = B H {t) + [ c{s)ds 
Jo 

is a fBm under the probability measure F H on F H defined by 

^^=exp^(-( W ,c))=cxp(- ^ c(s)dB H (s)--\\c\\%), 
where c is the the continuous function with supp c C [0, oo) satisfying 

c(s)(f)(s, t)ds = c(t) 



and 



Iff 



OO POO 



c(s)c(t)(f)(s, t)dsdt. 



to Jo 

Let us now multiply (fT4")) by the integrating factor 

P(t) :=cxp(- f b(s)ds), 
Jo 

so that we get 

d(j3(s)p(s) = P(s)a(s)ds + (3{s)q{s)dB H (s), 
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or by integrating 
(15) p(t)p(t) -- 

Assume that 

\\a\\ 2 ,: = Ep H 



I3(s)a(s)ds + J /3{s)q{s)dB H (s). 
a(s)a(t)<p(s, t)dsdt] 



OO fOO 



10 JO 

/* oo /*oo 

Ep H [ / D rt >(s)a(t)D' l '(t)a{s)dsdt} < oo, 
Jo Jo 



where denotes the </>-derivative with respect to B H . If we now apply the quasi- 
conditional expectation operator 



to ([15)) we get 
(16) 



p(t)v{t) = / p{a)E} H [a{a)\??]da 



From (|T6|) we get the solution 



(17) 



P(t) 



exp(- / &(r)dr)i? fJJ [a(s)|.Ff]ds. 



5. A INFINITE HORIZON MAXIMUM PRINCIPLE 



In this section we prove a maximum principle for systems driven by a fractional 
Brownian motion. For classical Brownian motion systems see e.g. |Hau86j . [PSOO], 
|Pen90| and |YZ99| and the refrences therein for more information. Assume we 
have an m-dimensional fractional Brownina motion, B (t), with Hurst parameter 
H = (Hi, H2, ■ ■ ■ , H m ). Let X(t) = X u (t) be a controlled fractional diffusion, 
described by the stochastic differential equation; 

dX{t) = b(t, X(t),u(t), u)dt + a(t, X(t),u(t), u)dB H ' (t) 

(18) X(0) = xeR n . 

Let 

£t c T u 

be a given subfiltration, representing the information available to the controller at 
time t;t > 0. The process u(t) is our control, assumed to be {£t}t>o adapted and 
with values in a set U C K™. Let A be our family of £ t -adapted controls. 
Let / : [0, 00] x l n x U x Q — >• R ra be adapted and assume that 



E 



\f(t,X(t),u(t),u)\dt 



< 00 for all u € 



Then we define 

J(u) = E 

to be our performance functional 



f(t,X(t),u(t),uj)dt 
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Definition 8 (Admissible pair). Let A denote the J-^-adapted processes u : 
[0,oo) x SI — > U such that X u (t) exists and doesn't explode in [0,oo) and such 
that it satisfies the stochastic differential equation (JTSJ). If u g A and X u (t) is the 
corresponding state process, we call (u,X u (t)) an admissible pair. 

We study the problem to find u 6 A such that 

J(u) = sup J(u). 

ueA 

If such u 6 A exists the u is called an optimal control and (u,X u (t)) is called an 
optimal pair. 

Let C([0, 00), W lXm ) be the set of continuous functions from [0, 00) into R" xm . We 
now define the Hamiltonian H : [0, 00) x R n x U x R" x C([0, 00), R nXm ) -> R, by 

i?(t, x, g, w) = /(t, a;, w, cj) + b T (t, x, u, uj)p 



n m „ c 
i=l i.k=l J ° 



(19) 



qi,k(s)())H k (s,t)ds, 



where <j>u k (s, t) is defined as above. For notational convenience we will in the rest 
of the paper suppress any u> in the notation. 

The adjoint equation that arrise in the maximum principle is the unknown Tt- 
predictable processes (p(t),q(t)) that satisfies the following stochastic differential 
equation; 

(20) dp(t) - -V x H(t,X(t),u(t),p(t),q(t))dt + q(t)dB H (t). 

For the solution of equations of the form (|17|) see Chapter 14.11 

Theorem 5.1 (Partial Information. Infinite Horizon Fractional Maximum Princi- 
ple). Let it £ A and let (p(t),q(t)) be an associated solution to the equation (|20p . 
Define 1 4 as 



i j,k 



OO POO 



Assume that for all u £ A the following terminal condition holds: 



(21) 
and 



< E 



lua[p(ty(X(t)-X(t))} 

t— too 



h < 0. 

Moreover, assume that 

H(t,x,u,p(t),q(t)) is concave in x and u and 



E 



(22) 

and 
(23) 



maxi? 

ueu 



H(t,X(t),u(t),p(t),q(t))\£ t 
H(t,X(t),u,p(t),Q(tWt 



E 



\V u H(t,X(t),u(t),p(t),q(t),f(t,-))f 



< 00. 



Then we have that u(t) is optimal. 
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Proof. Let I := E[f£°(f(t,X(t),u(t)) - f(t,X(t),u(t)))dt] = J(u) - J{u). Then 
/ = I± — I2 — I3 — I4, where 



h -=E 



(H(s,X(s),u( S ),p(s),q(s)) 



(24) -H(t,Jt(8)Mt)MM))d* 



p(sf '(6(a) - h(a))da 



/a := S 



X! J2 ( s ' u{s))% k {t)(j)H k (s, t)]dsdt 

i k 



where $i,fc(a, X(s), u(s), X(s), u(s)) := o^kpf (a), u(a)) - £7- i) &(X(s), w(s)). From 
concavity we get that 



(25) 



(26) 

Then we have from 



H(t, X(t), u(t), P (t),q(t)) - H(t, X(t),u(t), P (t),q(t)) 
< V x H(t, X(t), u(t),p(t),q(t)) T (X(t) - X(t)) 
+ V u H{t,X(t),u(t),p{t),q(t)) T (u(t) - u{t)). 



and that u(t) is adapted to £ , 



> V U E 



H(t,X(t),u,p(t),q(t))\£ t (u(t)-u(t)) 



u—u(t) 



(27) 



E V u J?(t, X(t), u(t),p{t), q{t)) T (u(t) - u(t))\£ t 



Combining ([20]) . ([24]) . ([26]) and ([27 



If < £ 

= £ 

We see that 

/a = £ 



V x H(t, X(s),u(s),p(s),q( S )) T (X(s) - X(s))ds 



(X(s) - X(s)) T dp(s) 



p(s) T (b(s) - b(s))ds 



-Ji. 





00 



£[/ p(*)(dX(t) - dX(t))] 



+ E[ / p T (i)K fc (X(a),w(a)) - a hk (X(s), u(s))}dB H (t)} 



E[ p(t)(dX(t) - dX(t))] 
Jo 
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Now, using ([2TTl and Ito's formula 



< E 

= E 

+ E 



lim[p(tY(X(t)~X(t))} 

t—^-oo 



(X(s) - X(s)) T dp(s) 



E[ p(t)(dX(t)-dX(t))} 



V V $i,fc(s, X(s),u{s),X(s),u(s))q^ k (t)<p Hk (s, t)]dsdt 

u 



% j,k 



i k 
poo poo 



Dp s {a iik (t,X(t),u(t)) - a ilk (t,X{t),u(t))}D*\q id {a)dtda] 



TOO 1 TOO I TOO 1 TOO 

— 1 2 ' J l ' J 3 ' J 4 ' 

Finally, combining the above we get 

J(u) - J(u) < 7f - / 2 °° - /3 00 - i| 

7-00 

'3 J 4 



^ TOO TOO TOO TOO 



< 0. 



This holds for all u £ A so the result follows. 



□ 
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We now give a infinite horizon version of the example in [BH0SO2J. 

Example 5.1 (minimal variance problem). Consider a financial market driven by 
two independent fractional Brownian motions, Bi(t) = B 1 1 (t) and B2(t) = B 2 2 (t) 
with Hi, H 2 & (5, 1), defined as follows: 

1: A Bond Price 

(dS (t) = 0, 
\S (0) = 1. 

2: Stock 1 

(dS 1 (t) = dB 1 (t), 

\S 1 (0) = 8 1 . 

3: Stock 2 

(dS 2 (t)=dB 1 (t)+dB 2 (t), 
\Sa(0) = «a. 

If4>{t) — (4>o(t), 0i (t) , 4>2{t)) € M 3 is a portfolio ( the number of units of bond, stock 
1 and stock 2, respectively, held at time t), then the corresponding value process is 

2 

V 9 (t) = 9(t)S(t)=J2^)Si(t)- 

j=0 

The protfolio is called self-financing if 

dV e (t) =6(t)dS(t) = 6i(t)dBx(t) + d 2 {t)(dB 1 (t) + dB 2 (t)). 

The market is called complete if any bounded J 711 -measurable random variable F, 
can be hedged, in the sense that there exist a self-financing portfolio 6{t) and an 
intial value z £ R such that 

/>oo 

F = z+ 9{t)dS(t), 
Jo 

almost surely (see [H0O3] for more details). Let us consider the case where we are 
unable to trade n stock 1. Let us say we want to stay "close" in some sense to Bi(t) 
at all times. We let 02(t) = u(t) and consider "close" as small L 2 (/i)- difference. 
Find z€i and admissible u(t) such that 

J(z,u):=E[{B 1 {T)-{z+ [ u(t)(dBi(t) + dB 2 {t))dt)} 2 ] 

Jo 

is minimal for all T > 0. Its clear that z = is optimal, so it remains to minimize 
J{z,u) :=E[{ ( (u(t) - l)(dBi(t) + dB 2 {t))dt} 2 ] 



Let 

J{u) = E 

where 



e-" l -X 2 {t)dt , 
dX(t) = u^dB^t) + (1 - u{t))dB 2 (t), 
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and p > 0. We have 

H(t,x,u,p,q) = \e~ pt x 2 



qx(s)4>x(s,t)ds + J q 2 (s)<j) 2 (s,t)ds 
q 1 {s)cj) 1 (s,t)ds. 

Jo 
So that 

V x H(t, x, u,p, q) = xe~ pt x 2 , 

and 

poo poo 

V u H(t,x,u,p,q) = q 1 (s)(j) 1 (s,t)ds + q 2 (s)<f> 2 (s,t)ds 
Jo Jo 

It now follows that we have 

dp{t) = e- pt X(t)dt + qx{t)dBx{t) + q 2 (t)dB 2 (t). 
If we let qi(t) — ^e^ pt u(t), q 2 (t) — ig~ pt (l — u(t)) then we have that 

H(t,X(t),v,p(t),q(t)) = \^ pt X 2 {t) 



+ 



V i.fo ~ erPS{l ^ ( t )1 ^'^ ds + J ^e- pt (l~u(t))s)(j> 2 (s,t)ds^ 



+ ( -e- ps u(s)Ms,t)ds. 
Jo P 

Since we need that the maximum of this expression is attained atv = u, is is a easy 
to see that we must have 

((1 - u(s))^i(s,i) - u(s)<p 2 (s,t))ds = 0. 

This is a symmetric Fredholm integral equation of the first kind and it is well known 
that it ha a unique solution u(t) 6 L 2 ([0, oo)), see [Fer09j. This u(t) satisfies all 
reqirements in our theorem and we also see from ltd 's lemma that 

p(t) = -e- pt X(t), 
P 

so that we have that 



lim p(t)(X(t) -X(t)) = 0, 

t— >oo 



almost surely. This gives us that 

u{t), 

is a optimal control. 
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